Path-wise solutions of SDEs driven by Levy 

processes 

David R. E. Williams* 
June 1999 



I 



Abstract 

In this paper we show that a path- wise solution to the following integral 
equation 

Yt= [ f(Yt)dXt Yo = ae^'' 
Jo 

exists under the assumption that Xt is a Levy process of finite p- variation 
for some p > 1 and that / is an a-Lipschitz function for some a > p. 
There are two types of solution, determined by the solution's behaviour 
at jump times of the process X, one we call geometric the other forward. 
The geometric solution is obtained by adding fictitious time and solving 
an associated integral equation. The forward solution is derived from the 
geometric solution by correcting the solution's jump behaviour. 

Levy processes, generally, have unbounded variation. So we must use a 
pathwise integral different from the Lebesgue-Stieltjes integral. When X 
has finite p- variation almost surely for p < 2 we use Young's integral. This 
is defined whenever / and g have finite p and g- variation for 1/p + l/g > 1 
(and they have no common discontinuities). When p > 2 we use the 
integral of Lyons. In order to use this integral we construct the Levy area 
of the Levy process and show that it has finite (p/2)-variation almost 
surely. 



Introduction 



In this paper I give a path-wise method for solving the foUowing integral equa- 
tion: 

Yt^Yo+ f f{Yt) dXt Yo = ae M'^. (1) 
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when the driving process is a Levy process. 

Typically, a Levy process a.s. has unbounded variation. The integral does 
not exist in a Lebesgue-Stieltjes sense. However, the integral still makes sense as 
a random variable due to the stochastic calculus of semi-martingales developed 
by the Strasbourg school 

The semi-martingale integration theory is not complete though. There are 
processes of interest which do not fit into the semi-martingale framework, for 
example the fractional Brownian motion. An alternative integral is provided by 
the path- wise approach studied by Lyons , jl^ and Dudley |^ . The basis of 
their papers is that of Young , who showed that the integral 

* / dg (2) 

is defined whenever / and g have finite p and g- variation for 1/p + 1/q > 1 
(and they have no common discontinuities). For a comprehensive overview of 
the theory we recommend the lecture notes of Dudley and Norvaisa B . 



Recently in 1 13 , a system of linear Riemann-Stieltjes integral equations is 



solved when the integrator has finite p- variation for some < p < 2. These 



results are contained in Theorem 1.1 where we allow non-linearity of the vector 
field /. This is because our approach is an extension of the method of ||], JTof . 

The approach that I follow distinguishes two cases. The first is when the 
process has finite p- variation a.s., for some p < 2. We use the Young integral 
[ p^ . In 1^ is solved when Xt is a continuous path of finite p- variation for 
some p < 2. 

The second case is when the process has finite p-variation a.s., for some 
p > 2. The Young integral is only defined when / and g have finite p and 
g- variation for l/p+ l/q > 1. So an iteration scheme on the space of paths with 
finite p-variation does not work. However, Lyons defined an integral against a 
continuous function of p- variation for some p > 2 The integral is developed 
in the space of geometric multiplicative functionals (described in Appendix]^. 
The key idea is that we enhance the path by adding an area function to it. If 
there is sufficient control of the pair, path and area, then the integral is defined. 
The canonical example in |l0[| is Brownian motion. The area process enhancing 
the Brownian motion is the Levy area Ch.7, Sect. 55]. I show that there is 
an area process of a Levy process which has finite (p/2) -variation a.s.. 

In order to solve (|l|) for a discontinuous function I add fictitious time during 
which linear segments remove the discontinuities, creating a continuous path. 
By solving for the continuous path and then removing the fictitious time we 
recover a solution for the discontinuous path. This is called a geometric solution. 
A second type of solution is derined from the geometric solution which we call 
the forward solution. 

The first section treats the case where the discontinuous driving path has 
finite p- variation for some p < 2. The second section treats the case where the 
path has finite p- variation for some p > 2 only. The main proofs of the second 
section are deferred to the third section. In the appendix I prove the homeo- 
morphic fiow property for the solutions when the driving path is continuous. 
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This is used in proving that forward solutions can be recovered from geometric 
solutions. 

1 Discontinuous processes - p < 2 

In this section we extend the results of |^ to allow the driving path of (|^) to have 
discontinuities. The results are applied to sample paths of some Levy processes, 
those that have finite p-variation a.s. for some p < 2. Throughout this section 
p & [1,2) unless otherwise stated. 

First, we determine the solution's behaviour when the integrator jumps. 
There are two possibilities to consider: the first is an extension of the Lebesgue- 
Stieltjes integral; the second is based on a geometric approach. 

Suppose that the discontinuous integrator has bounded variation. The solu- 
tion y would jump 



at a jump time t of x. If x has finite p- variation for some 1 < p < 2 we insert 
these jumps at the discontinuities of x. We call a path y with the above jump 
behaviour a forward solution. 

The other jump behaviour we consider is the following: When a jump of 
the integrator occurs we insert some fictitious time during which the jump is 
traversed by a linear segment, creating a continuous path on an extended time 
frame. Then we solve the differential equation driven by the continuous path. 
Finally we remove the fictitious time component of the solution path. We call 
this a geometric solution because the solution has an 'instantaneous flow' along 
an integral curve at the jump times. This jump behaviour has been considered 
before by |ll| and [|. 

The disadvantage of the first approach is that the solution does not, gener- 
ally, generate a flow of diffeomorphisms 0. 

In this section we prove the following theorem: 

Theorem 1.1 Let Xt he a discontinuous function of finite p-variation for some 
p < 2. Let f be an a-Lipschitz vector field for some a > p. Then there exists a 
unique geometric solution to the integral equation 



With the above assumptions, there exists a unique forward solution as well. 

Before proving the theorem we recall the definitions of p-variation and a- 
Lipschitz: 

Definition 1.1 The p-variation of a function x{s) over the interval [0,t] is 
defined as follows: 



yt - yt- = f{yt-) [xt - Xt-) 





(3) 




sup 

7re7r[0,t] 
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where 7r[0,i] is the collection of all finite partitions of the interval [0,t]. 

Remark 1.1 This is the strong p-variation. Usually probabilists use the weaker 
form where the supremum is over partitions restricted by a mesh size which tends 
to zero. 

Definition 1.2 A function f is in Lip(ck) for some a > 1 if 

df 

ll/lloo < oo and — S Lip{a - 1) j ^l,...,d. 



Lip(a) = ll/IU+Ell|rllLip(a-l) for a > 1. 



Its norm is given by 

d 

I Qrj. llLip(Q-l) 

j=l 3 

This is Stein's definition of a-Lipscliitz continuity for a > 1. It extends 
tlie classical definition: / is in Lip(Q;) for some a € (0, 1] if 

\f{x)-f{y)\<K\x-yr 

with norm 

1.1 Geometric Solutions. 

In this subsection we define a parametrisation for a cadlag path x of finite p- 
variation. The parametrisation adds fictitious time allowing the traversal of the 
discontinuities of the path x. We prove that the resulting continuous path x^ 
has the same p-variation that x has. We solve driven by using the method 
of Lyons §. Then we get a geometric solution of (||) by removing the fictitious 
time (i.e. by undoing the parametrisation). 

Definition 1.3 Let x be a cadlag path of finite p-variation. Let (5 > , for 
each n>l, let tn be the time of the n'th largest jump of x. We define a map 
: [0,T] [0,r + ^X^i^i (where j{u) denotes the jump of the path x 

at time u) in the following way: 

oo 



T\t)=t + 5^\j{t^)\PX{u.<t}{t). (4) 



n=l 



The map : [0,T] — > [0,t^(T)] extends the time interval into one on where we 
define the continuous process x^ {s) 

xHs) 

_( Xit) ifs^T^it), 

- 1 x{t-) + [s THt-))j{t^)5-^\j{tn)\-^ ifs e [T'{t-y{t,,)). 
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X(s) / 


> 










\ ^^^^ 

^ > 





T(s) 



Figure 1: Let Xs be a discontinuous path of bounded variation {p = 1). De- 
fine a map r(s) inserting fictitious time for tlie discontinuities of x. Define a 
parametrisation (i(s), j/(s)) in the manner of (|^). {t{s),y{s)) traverses the jumps 
of X during the fictitious time. 



Remarks 1.1 

1. (s, xf), s e [0, r*(r)] is a parametrisation of the driving path x. 

2. The terms \j{tn)\^ m (|) ensure that the addition of the fictitious time 
does not make T^{t) explode. 

3. In Figure |^ we see an example of a parametrisation of a discontinuous 
path Xs in terms of the pair {t{s), y{s)). 

The next proposition shows that the above parametrisation has the same p- 
variation as the original path, on the extended time frame [0, {T)]. 

Proposition 1.1 Let x be a cddldg path of finite p-variation. Let x^ be a 
parametrisation of x as above. Then 

l|a;'^llp,[o,r-5(T)] = l|a;||p,[o,T] > 0. 

Proof. Let ttq be a partition of [0,T''(r)]. Let 
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We show that we increase the value of Vp{'!To) by moving points lying on the 
jump segments to the endpoints of those segments. 

Let ti^i,ti,tijf.i be three neighbouring points in the partition ttq such that 
ti lies in a jump segment. Consider the following term: 

+ -<r (6) 

We show that is dominated by replacing xf, by one of and x^, where I 
and r denote the left and right endpoint of the jump segment containing ti. 
For simplicity we set a = , b — xf,^_^ and c = xf . Let 

L = {c + kx : k e {0,1), c,x eK.'^,xj^ 0}, a,be m.''\L. 

Let the function / : [0, 1] (0, oo) be defined by 

f{k)^\a-d\P + \d-b\P ,d^c + kx. 

Then / e C'^[0, 1] and one can show that /" > on (0, 1) when p > 1. To 
conclude the proof we move along the partition replacing ti which lie in the jump 
segments by new points t'i that increase Vxs{ttq). The partition ttq is replaced 
by a partition ttq whose points lie on the pre-image of [0,r''(r)]. Therefore we 
have 

Va:s{TTo) < V^s{tt'q) = 14(7ro). 
Hence ||a;''||p,[o,r'5(T)] = lk||p,[o,T]- □ 

Theorem 1.2 Let x be a cadlag path with finite p-variation for some p < 2. 
Let f be a Lip(7) vector field on R" for some j > p. Then there exists a unique 
geometric solution y, having finite p-variation which solves the differential equa- 
tion 

dyt^f{Vt)dxt 2/0 = a eR". (7) 

Proof. Let x^ be the parametrisation given in (||). The theorem of section 
three of ||^ proves that there is a continuous solution which solves on 
[0, T^{T)\. Then (s, ) is a parametrisation of a cadlag path y on [0, T]. 

The solution is well-defined. To see this, consider two parametrisations of x 
and note that there exists a monotonically increasing function As such that 

{S,xi)^{\s,xl^). □ 

1.2 Forward Solutions. 

In this subsection we show how to recover forward solutions from geometric 
solutions. The idea behind our approach is to correct the jump behaviour of the 
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geometric solution using a Taylor series expansion Lemma |l.l| . The correction 
terms are controlled by 



oo 



J2 1^*- " ^t- P 

1=1 

which is finite due to the finite p- variation of the path x. 

In the case where the driving path has only a finite number of jumps we note 
that the forward solution can be recovered trivially. It is enough to mark the 
jump times of x and solve the differential equation on the components where x 
is continuous, inserting the forward jump behaviour when the jumps occur. It 
remains to show that the forward solution exists when the driving path has a 
countably infinite number of jumps. The method we use requires the following 
property of the geometric solution: 

Theorem 1.3 Let x be a continuous path of finite p-variation for some p > I. 
Let f be in Lip(a) for some a > p. The maps {'!Tt)t>o ■ K" obtained 

by varying the initial condition of the following differential equation generate a 
flow of homeomorphisms: 

d-Kt — /(ti"*) dxt ttq = /rf, {the identity map). (8) 



We leave the proof of Theorem 1.3 until Appendix |^. We note the uniform 
estimate 

sup < C{T)\a-b\. (9) 

0<t<T 

The following lemma will enable estimates to be made when the geometric 
jumps are replaced by the forward jumps: 

Lemma 1.1 Let x be a cddldg path with finite p-variation. Let f be in Lip(a) 
for some a > p. Let Ayi (resp. lS.Zi) denote the geometric (resp. forward) 
solution's jump which correspond to Ax^, the i'th largest jump of x. Then we 
have the following estimate on the difference of the two jumps: 

||Ay, - A2,||oo <K\Ax^\'^ 
where the constant K depends on ||/||Lip(Q)- 

Proof. Parametrise the path x so that it traverses its discontinuity in unit 
time. Solve geometrically over this interval with the solution having initial 
point a. Note that the forward jump is the first order Taylor approximation to 
the geometric jump. Then 



yi{aj^yo{a) + 



ds 
zi{a) - 



1 d^y^ia) 



.=0 2 ds^ 

1 d'^ysja) 

2 ds^ 



for some < 6 < 1 

(10) 
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We estimate the second order term by 



1 d'^Vsia) 




2 ds2 


OO 



1 d_ 

2 lis 



/(y.(a))(Ax,) 



< - ||V/|| 11/11 \Ax 



|2 

iLip(a) I 



< ^\\f\\:.___,_^^x,'^ 



(11) 



Both ||V/||oo and ||/||oo are finite because / is Lip(a) for some a > p > 1. 
□ 

Theorem 1.4 Let x be a cadldg path with finite p-variation. Let f be in Lip(a) 
for some a > p. Then there exists a unique forward solution to the following 
differential equation: 



dzt = f{zt) dxt 



zq = a. 



(12) 



Proof. By Theorem 1.3 there exists a unique homeomorphism y which solves 
dyt = f{yt) dxt 2/0 = a 

in a geometric sense. 

Label the jumps of x by jx = {jiji^i according to their decreasing size. Let 
z" denote the path made by replacing the geometric jumps of y corresponding 
to {ji}?=i by the forward jumps {/() (Axi)}"^^. We show that the {z")n>i have 
a uniform limit. 

We order the corrected jumps chr onologically, say Then we estimate 

the following term using Lemma 1.1 and the uniform bound on the growth of y 
given in (||): 

n 

W:{a) ys{a)\ < - yuAyu-uuK-A'^M 

n 

< ciT)Y,Wui<^)-yu-,MK-A^))\ 

oo 

< C\T)KJ2\^^^\^■ (13) 

i=l 

So we have the uniform estimate 

oo 

lk"-2/||oo<if(C3(T),||/||Lip(a))^|Aa;,|2<cx) Vn > 1. (14) 

1=1 

We use an analogous bound to get Cauchy convergence of {z"}n>i- Let 
m, r > 1. 



<i^(C(T,z™),||/||Lip(„)) £ I A. 



i—7n~\'l 
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One notes that {C(T, z"^)} are uniformly bounded, because of the boundedness 
of C(T) = C{T,y) and the Lipschitz condition on /. Therefore we have the 
following estimate: 

oo 

This implies that {z"} are Cauchy in the supremum norm because x has finite p- 
variation {p < 2) which implies that J2m+i l^^;^^ tends to zero as m increases. 
□ 



Remark 1.2 Theorems and LJ combine to prove Theorem 



1.1. 



Corollary 1.1 With the above notation, z has finite p-variation. 
Proof. Let s < < e [0,T]. 

\zt - Zs\<\{zt- Zs) - [yt -ys)\ + \yt - ys\ 

where (yt — ys) is the increment of the geometric solution starting from z^ driven 
by the path xt on the interval [s,T]. Then 

\{zt - Zs) - {yt - ys)\ < C ^ |Axjp and \yt - ys\ < \\x\\p^[s,t] , 



which implies that 



\zt-Zs\^<2P-'i^C''{ Y: \Ax.\T + MUs,U' 



Jx|[s,f] 



hence 



l^llp,[0,T] 



1/p 

< oo. 



Jx |[0,T] 

□ 



1.3 p-variation of Levy processes 



In this subsection we apply Theorem 1.1 to Levy processes which have finite 
p-variation a.s.. 

Levy processes are the class of processes with stationary, independent incre- 
ments which are continuous in probability. The class includes Brownian motion, 
although this process is atypical due to its continuous sample paths. Typically 
a Levy process will be a combination of a deterministic drift, a Gaussian process 
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and a jump process. For further information on Levy processes we direct the 
reader to |l| . 

The regularity of the sample paths of a Levy process has been studied in- 
tensively. In the 1960's several people worked on the question of characterising 
the sample path p-variation. The following theorem, due to Monroe, gives the 
characterisation: 

Theorem 1.5 J7^ , Theorem 2] Let {Xt)t>o he a Levy process in R" without a 
Gaussian part. Let v he the Levy measure. Let (3 denote the index of Xt, that 
is 

/3 = inf|a>0:^ \y iy{dy) < ooj (15) 
and suppose that 1 < f3 < 2. If ^ > [3 then 

^{\\X\l,<^^=l (16) 
where the ^-variation is considered over any compact interval. 



Remark 1.3 Note that all Levy processes with a Gaussian part only have finite 
p-variation for p > 2. 



Corollary 1.2 Let {Xt)t>o he a Levy process with index (3 < 2 and no Gaussian 
part. Let f he a vector field in Lip(a) for some a > p. Then, a.s., the following 
stochastic differential equation has a unique forward and a unique geometric 
solution: 

dYt = f{Yt) dXt Yo = a. 



Proof. The corollary follows immediately from Theorems 1.5 and 1.1 



□ 



2 Discontinuous processes - p > 2 



The goal of this section is to extend (Corollary 1.2) to let any Levy process be 
the integrator of (|l|). 

One problem we have is that the Young integral is no longer useful because 
we use a Picard iteration scheme which fails condition (0) when p > 2. However, 



we can use the method from |10 . To define the integral we need to provide more 
information about the sample path. We do this by defining an area process of 
the Levy process. Then we prove that the enhanced process (path and area) 
has finite p-variation Definition |A.3| . 

We parametrise the enhanced process in an analogous manner to (|) (adding 
fictitious time). Then we solve (|) in a geometric sense using the method for 
continuous paths {p > 2) given in [ p^ . Finally, forward solutions are obtained 
by jump correction as before. 
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Before enhancing {Xt)t>o we give an example which shows that there exist 
Levy measures with index two. So a Levy process does not need a Gaussian 
part to have, a.s., finite p- variation only for p > 2. 

Example 2.1 One can define the following measures onR: 

Vk {dx) = Ixp^+i/fe dx \x\ e ((fc + l)-3('=+i),/fc-3fe] A 

r]m (dx) = ^ i^k {dx n Jfc n (- Jfc))- 
fc=l 

We show that rj = limm^oo "Hm is a Levy measure. The integrability condition 

\x\^-q {dx) < oo. (17) 



must be satisfied. 

\x\'^r]m {dx) 

\x\<l 



pi m m 

•^0 k=l k=l 
2f]fc|fc-3_(fc + i)-3(i+iA)j 

k=l ^ ' 

2 ^ fc J _ 2-3(l+l/fc) ^-3(l + l/fc) 

k=l 

= 2f]fc-2|i_2-3(i+i/'^-)fc-3/4 

k=l ^ ^ 

oo 

CXD, 



kx^/"" 



(fc+l)-3(fc+l) 



< 



fc=l 



where C is some suitable constant. We take the limit as m tends to infinity on 
the left hand side to prove (pi 
Now we show that 



\x\°'r] [dx) — oo 

\x\<l 

for all a < 2. Fix a <2. Define the following number: 

m{a) = infjfc : a + 1/k < 2} < oo as a < 2. 
Let m > m{a). Then 



(18) 



\x\"rirn {dx) 



\x]<l 
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> 2 1 /^^-3fc(a + l/fc-2) _ , i\-3(k+l)(a+l/k-2)\ 

> ^ ^ |(^^ _^ -^^3(fc+l)(2-(a+l/fe)) _ ^3fc(2-(a+l/fc))| 



as m — > cxD . 



rftjs proves that the index P of rj equals two. Theorem L5 implies that the pure 
jump process associated to the Levy measure rj a.s. has finite p-variation for 
p > 2 only. 

The following theorem gives a construction of the Levy area of the Levy 
process {Xt)t>o. The Levy area process and the Levy process form the enhanced 
process which we need in order to use the method of Lyons [ p^ . 

Theorem 2.1 The d-dimensional Levy process {Xt)t>o has an anti- symmetric 
area process 

{As,tr ^\j\l_odXi-Xi_odXl z,j = l,2. a.s. 
The proof is deferred to Section ^. 

Theorem 2.2 The Levy area of the Levy process {Xt)t>o a.s. has finite (p/2)- 
variation for p > 2. That is 

sup(5:|A,_.,r/^)^/^ <^ a.s. 

7T 

where the supremum is taken over all finite partitions tt of [0, T] . 

The proof is deferred to Section ||. 

Now we parametrise the sample paths of {Xt)t>o as before (||). 

Proposition 2.1 Parametrising the process {Xt)t>a does not affect the area 
process' (p/2) -variation. 



Proof. The proof is similar to the proof of Proposition 1.1. One can show that 
if A lies in a jump segment then 

is maximised when A is moved to one of the endpoints of the jump segment. 
□ 

With the parametrisation of the path and the area we can define the integral 
in the sense of Lyons . Consequently we have the following theorem: 
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Theorem 2.3 Let {Xt)t>o be a Levy process with finite p-variation for some 
p > 2. Let f be in Lip(Q!) for some a > p. Then there exists, with probability 
one, a unique geometric and a unique forward solution to the following integral 
equation: 



Yt = Yo+ f f{Yt) dXt Yo^ae 
Jo 



(19) 



Remark 2.1 When constructing the forward solution it is necessary that the 
sum 

oo 

|2 



n=l 

remains finite. This is guaranteed by the requirement on Levy measures to satisfy 

|xp A 1 iy{dx) < oo. 

x\<l 



3 Proofs of Theorem 12. 11 and Theorem 12.2 



For clarity throughout this section we assume that the Levy process {Xt)t>o is 
two dimensional and takes the following form: 

Xt = Bt+ [ X (Ntidx) - tv{dx)). (20) 

J\x\<l 

That is, {Xt)t>() is a Gaussian process with a compensated pure jump process, 
whose Levy measure is supported on (a; e : |a;| < 1). 

Proposition 3.1 The d-dimensional Levy process {Xt)t>o has an anti- symmetric 
area process 

(As.tY' =IJ' ° dXi - Xi^ o dXl z, J = 1, 2. a. s. 
For fixed s < t we obtain the area process by the following limiting procedure: 

n 2'"-l 



n — *oo 

m=0 fe=l, 
odd 



where A^j!^ is the area of the (ij) -projected triangle with vertices 

X{U(^k+l)/2,m~l), ^(U(fc_l)/2,m-l)j X{uk,m) 

where Uk^m — s + fc2^™(f — s). Also we have the second order moment estimate 

{A'Lr)<C{iyKt-sr. (21) 
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Proof. We define Ag^tin) 

1 2"-l 

2"-l 
fc=0 

where -Bfe^„ is tlie (signed) area of the triangle with vertices 

X{s), X{uk,n), X{uk+l,n)- 

By considering the difference between Asotin) and Asotin + 1) we see that 

B2k,n+1 + B2k+l,n+l — Bk,n 

is the area of the triangle with vertices 

X{uk,n), X{uk+l,n), -'^(M2fe+l,n+l) 

which we denote by Ak,n- We re-order As^t{n) 

n 2™-! 

2 



m=0 fe=l, 
odd 



2 " ' " / 

® ^^(W(fe+l)/2,m-l) - -'^(W(fe-l)/2,m-l)^ 



^ n 2"-! 

2 E E 



m=0 fe=l, 

where dfe,^ = 1/2 {X{u(^k+i)/2,m-i) + X{u(^k-i)/2,m-i))- The convergence to 
the area process is completed using martingale methods. 
Let dn = (r{X{uk,n) ■■ A; = 0, . . . , 2") . Then 

Lemma 3.1 

K(^X{uk,m) \ dm-i^ = dk,m a. s. (22) 
Proof. For ease of presentation we let 

Ui = X{Uk,m) - ^(W(fe-l)/2,m-l) 
U2 = X{{U(^k+l)/2,m-l) - X{Uk,m)- 
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Then 

,m I — 1 

= E^X{uk^m) - dk,m I -'^('«(fc-l)/2,m-l), -^("(fc+l)/2,m-l)^ 
= 2 ^(Ul ~ U2 I X{u(^k-l)/2,m-l),X{u(k+l)/2,m-l)^ 

Using the stationarity and the independence of the increments of X we see 
that Ui and U2 are exchangeable, that is 



Ui eA,U2e Bj = F\^U2e A, Ui e Bj VA, B e <B(]R^). 
The exchangeabihty extends to the random variables 

{Ui\X{u(i,-l)/2.7n-l), ^(W(fe+l)/2,m-l)) « = 1, 2. 

We deduce that 

E^C/i - U2 I X{u(^k-l)/2,m-l),X{u(^k+l)/2,m-l)^ =0. □ 



Returning to the proof of Proposition 3.1 , we compute the variance of A^^; 
This will be used to show that 

sup E( As t{nf ) < oo. 
n>l V ' / 



= e(([x(^)K,„) - 4i][{/f ) + c/f ] - - d^l][u\^^ + 

^ (1) + (2) + (3). 

We use the independence of the increments and Ito's formula for discontinuous 
semi-martingales to compute (1), (2) and (3). 
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By applying Ito's formula and using the stationarity of the Levy process we find 
that 

(3) = (1) = 2-2™ {t - sf [ u{dx) [ Ixal' u{dx). 

J\x\<l J\x\<l 

Another application of Ito's formula gives 

(2) = -2E(^C/f ^C/f f/(^)c/f = -2E(ui'^u['^y(u^'^ui''' 
/ r \ 2 

-2m+l /J. \2 



= -2-^™+^ {t - s)M / xi X2 v{dx) . 

\J\x\<l J 

Collecting the terms together we have the following expression: 



where 



Co{y) = { / \xi\^ v{dx) I \x2\^ ^{dx) - [ I xiX2i'{dx) 

kJ\x\<1 J\x\<1 V^|x|<1 

Now we estimate the following term: 

n 2"-! s 2 



E(Al,in))=E((Y: E 

^ ^ ^^m=l fe=l, 



odd 

which through conditioning and independence arguments equals 

n 2"'-l 



m=l k=l, 
odd 



m=l fc=l, 
odd 

oo 2"^~1 



< c^oM E E - = c{t^){t - sf 



m=l fc=l, 
odd 

Wc use the martingale convergence theorem to dediice that a.s. there is a unique 
limit of Ag^tin). Furthermore the last calculation implies that there is a moment 
estimate of the area process given by 

kIa') <C{u){t-sf. □ 
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We note that there is another way that one could define an area process of a 
Levy process. One could define the area process for the truncated Levy processes 
and look for a limit as the small (compensated) jumps are put in. Using the 
above construction one can define Al ^ for a fixed pair of times, corresponding 
to the Levy process X^. With the cr-fields (25'^)e>o defined by 

6^ ^ (j{X^ : (5 > e) for e > 

we have the following proposition: 

Proposition 3.2 {A'^^^)^yQ form a {(&'^) -martingale. 

Proof. Let 77 > e > 0. By considering the construction of the area given above 
for the truncated processes X'' and X'^ we look at the difference at the level of 
the triangles A^ ^ and A^ ^ . 

where the superscript 77, e signifies that the process is generated by the part of 
the Levy measure whose support is (e, 77]. Using the spatial independence of the 
underlying Levy process we have 

= e(a^-) +e((^^::-C)) ®(^oc+i)An-i-^(i-i)/2,«-i) 
= 0. □ 

With the uniform control on the second moment of the martingale 
e({AIaA <C{^){t~sr Ve>0 



we conclude that Al ^ converges a.s. as e ^ 0. 
The algebraic identity 

As,u = + + ,J s<t<u (23) 

for the anti-symmetric area process A generated by a piecewise smooth path X 
extends to the area process of the Levy process. This is due to ( |2^ ) holding for 
the area processes A"^ of the truncated Levy processes X"^. 
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Proposition 3.3 The Levy area of the Levy process {Xt)t>o has finite {p/2)- 
variation for p > 2 a.s.. That is 



|p/2\2/p 



< oo 



a. s. 



where the supremum is taken over all finite partitions tt of [0,T]. 



Proof. In Proposition 3.1 we constructed the area process for a pair of times, 
a.s.. This can be extended to a countable coUection of pairs of times, a.s.. In 
the proof below we assume that the area process has been defined for the times 



fc2-"r, (fc + i)2-"r 



0,1,. 



,2"-l, 



n> 1. 



The proof follows the method of estimation used in To estimate the 
area process for two arbitrary times u < v we split up the interval [u, v] in the 
following manner: 

We select the largest dyadic interval [{k — 1)2^"T, fc2^"] which is contained 
within Then we add dyadic intervals to either side of the initial inter- 

val, which are chosen maximally with respect to inclusion in the interval 
Continuing in this fashion we label the partition according to the lengths of the 
dyadics. We note that there are at most two dyadics of the same length in the 
partition which we label [li^k,'fi,k] and [l2,k,f2,k] where ri^k < l2,k- Then 

oo 

[u,v] = IJ IJ [h,k,n,k]- 

fc=l i=l,2 

We estimate Au,v using the algebraic formula (E^). 



A, 



fc=l i=l,2 

l<a<6<2 l<j<k<m 



Noting that 



E E I 

l<a<fc<2 l<j<fc<m 
= E E \iXr...-Xl.^.)®iXr 



l<a<b<2 l<j<k<m 

^ E E 

l<a<6<2 l<j<k<m 



^ k=l i=l,2 
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we have the estimate: 



fc=l 4=1,2 



k=l i=l,2 



(24) 



Using Holder's inequaUty, with p > 2 and 7 > p — 1, we have 

(p/2)-l 00 . .p/2 



^„-7/((p/2)-l) 



n=l ^1=1,2 ^ 

00 

<ci(p,7) E"' E 

00 

+ C2(p,7) E"' E 



p—i 00 y 

71=1 ^i = l,2 



(25) 



ri=l i=l,2 



One can uniformly bound for any pair of times u < v Cz [0,T] by 

extending the estimate in ( p5| ) over all the dyadic intervals at each level n, that 
is, 

00 2" 

K..r/'<ci(p,7) E"'Ei^'.-'--i'^' 

n—l i—1 

00 2'^ 

+ C2(p,7) E""Ei^'-.'=-^'.J'- 

If the right hand side is finite a.s. then the area can be defined for any pair of 
times. 

The (p/2)-variation of the Levy area can be estimated by the same bound. 

00 2" 

sup^i^„,.r/2 < ci(p,7) E«"Ei^'.-'-.^-i'^' 



n—l i—1 



+ C2(p,7) E""Ei^'-.'=-^'.^-i'- 

n=l 2=1 

We use ( ^l| ) to control the first sum 

E(^|A,,t|P/2^ <C (t- s)P/2 for p < 4. 



(26) 



So we have 



E E E i^'.-'-.^ 1'^' < ^ E E(2""?^)"^ 



n—l i— 1 



n—l 2—1 
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OO 
Tl=l 

< cxD for p > 2. 

This imphes that the first term in the right hand side of (^6|) is a.s. finite. Now 
we consider the second term of (|2^). 

Lemma 3.2 

OO 2"-l 
n=l k=0 

Before proving the lemma we recall a result of Monroe . 

Definition 3.1 Let Bt be a Brownian motion defined on a probability space 
A stopping time T is said to be minimal if for any stopping time 

S < T, B{T) B{S) implies that a.s. S = T. 

Theorem 3.1 JT^ , Theorem 11] Let {Mt)t>o be a right continuous martingale. 
Then there is a Brownian motion (fi, 0f,_Bf) and a family (Tt) of &t- stopping 
times such that the process BTt has the same finite distributions as Mf. The 
family Tt is right continuous, increasing, and for each t, Tt is minimal. More- 
over, if Mt has stationary independent increments then so does Tt. 



Remark 3.1 It should be noted that the stopping times Tt are not generally 
independent of Bt. However, in the case of a- stable processes < a < 2 one 
can use subordination to gain independence of the stopping times 



Proof of Lemma |3.2| Let {Tt)t>o denote the collection of minimal stopping 
times for which 

Xt =^ Brt . 

The proof will be completed once it has been shown that 

OO 2"-l 

^n'' ^ |B^((fc+i)2-"T) - ■Br(fc2-"T)r < OO a. s. (27) 

n=l k=0 

The following inequality holds because Brownian motion is (l/p')-H61der con- 
tinuous a.s. for p' > 2: 

|5r((t.+,,„) -Sr(U.„)r < C|r(ife+i,„)-r(tfe,„)|^ (28) 

Vfc = 0,. .. ,2" - l,Vn > 1, a.s. 

where tk,n = fc2-"r and 2 < p' < p. 

15, Theorem 1] shows that the index of the process t(s) is half that of the 
Levy process. Therefore, with probability one, t(s) has finite (1 + (5)-variation 
for all 5 > 0. 



Path-wise solutions and Levy processes 



21 



Theorem 3.2 \I2, Theorem 5] If t is a minimal stopping time andK^BT-) = 0, 
then E(t) = E{B^. 

Consequently the process {Tt)t>Q can be controlled in the following way: 



Ttj^Ei^B',^j=E\^Xt'j=tJ^^ \x\'iy{dx) (29) 

where f is the Levy measure corresponding to the process Xt- From ( p9| ) and 
Theorem 3.1 we note that the process Tf is a Levy process whose Levy measure, 
say /i, satisfies the following: 

X ^{dx) < oo. 



From this result we deduce that the process Tt a.s. has bounded variation. From 
[|6[ Theorem 5] we note that there is a positive constant A such that 

Tt<At,yt>0^=l. 

From the above bound and using the fact that r has stationary independent 
increments one can show 

T{tk+l,n) - T{tk,n) < A{tk+l.n " tk,n) = ^2"" | r(tfc,„)j = 1, 

2"-l 



= 1. 



Returning to (Eq) we see that 



|i3r((t,+i,„) - Br(t.,„)r < C |r(tfe+i,„) - Titk,n)\ 



< CA2-"(f) 



which implies that 

oo 2"-l 



7o-"(4-i) 



Y.n-'Y. \Brak+i)2-^T) - S.(fe2-"T)r <CAJ2 nn-''^V-^^ < oo 



due to p' being chosen in the interval (2,p). □ 

This lemma concludes the proof that the bound in (^6|) is finite, which shows 
that the area process a.s. has finite (p/2)-variation. □ 

In this section we have proved that the area process exists and has finite 
(p/2)-variation when {Xt)t>a has the form ([20|). To prove Theorems 2.1, 2.2 we 
note that a general Levy process has the form 



Xt = at + Bt + Lt+ 



a. s. 



0<s<t 
\AX^\>1 
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So, we need to add area corresponding to the drift vector and the jumps of size 
greater than one. However, this part of the Levy process has bounded variation 
and is piecewise smooth so there is no problem defining its area. Similarly, it 
has a.s. finite (p/2)-variation. 

A Homeomorphic flows 

In this section we give a proof that the solutions, generated by (||) as the initial 
condition is varied, form a flow of homeomorphisms when the integrator is a 
continuous function. The proof modifies the one given in for the existence 
and uniqueness of solution to (Q). The main idea is that one uniformly bounds 
a sequence of iterated maps which have projections giving the convergence of 
the solutions with two different initial points and bounding the difference of the 
solutions. 

First, we need some notation. 

Definition A.l Let r("-'(R'^) denote the truncated tensor algebra of length n 
over M.'^ . That is 

n 

i=0 

where (M'^)®° = M and r(°°)(R'*) denotes the tensor algebra over R''. 

Let A = [0, T] X [0, T]. A map X : A ^ T(")(M'') will be called a multiplica- 
tive functional of size n if for all times s < t < u in \Q,T] the following relation 
holds in T(")(R'') : 

Xst ® Xtu = Xsu 

andxit'' = 1. 

A map X : A r*^")(R'') is called a classical multiplicative functional if 
t ^ Xt ^ x'^^ is continuous and piecewise smooth and 

X^ ^11 dX^, . . . dX^^ (30) 

J J S<Ul<...<Ui<t 

where the right hand side is a Lebesgue-Stieltjes integral. We denote the set of 
all classical multiplicative functionals in T(")(R'') by S'''"'(R''). 

Definition A. 2 We call a continuous function w : A — > R+ a control function 
if it is super-additive and regular, that is, 



uj {s,t) + Lu {t,u) < uj {s,u) V s < t < -u e [0, T] 
uj{s,s) = Vse[0,T]. 
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Example A.l Let X be a path of strong finite p-variation. Then we can define 
the following control function: 

= (31) 



Definition A.3 A functional X = {1,X^^\... defined on r(")(R'') 

where n — [p] is said to have finite p-variation if there is a control function 
CO such that 

|X«|<^^^^ V(,s,i)GA,z = l,...,n (32) 
for some sufficiently large [3 and x\ = T{x + 1). 

Theorem A.l J7^ , Theorem 2.2.1] Let X*^"^ be a multiplicative functional of 
degree n which has finite p-variation, with n ^ [p] ([p] denotes the integer part of 
p). Then for m > n there is a unique multiplicative extension X^™-* in T^"^\W^) 
which has finite p-variation. 



Remark A.l The above theorem shows that once a sufficient number of low 
order integrals associated to a path Xt have been defined, then the remaining 
iterated integrals of Xt are defined. 



Definition A. 4 We call a multiplicative functional X : A ^ T'^")(M'^) geomet- 
ric if there is a control function lu such that for any positive e there exists a 
classical multiplicative functional Y{e) which approximates X in the following 
way: 



Ystie) 



<eLjis,ty/P i = l, 



We denote the class of geometric multiplicative functionals with finite p 
-variation by f}G(IR'')P. 



Example A. 2 Let Wt be an M. -valued Brownian motion. Then the following 
functional W defined on T(2)(]R'i) belongs to nG{R'^)P for any p > 2. 

Wst = (l, Wt - Ws, II odWu, o dWu^ (33) 

V J J S<Ui<U2<t J 

where odWu denotes the Stratonovich integral. It should be noted that if one 
replaced the Stratonovich differential in (|3^) by the ltd differential then one 
would not get an element of ilGiW^Y . This is due to the quadratic variation 
term which occurs in the symmetric part of the area process 

^11 dWu,dWu2. 
J J S<Ui <U2 <t 
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It was shown in that one had sufficient control of the above functional 
to generate path-wise solutions to SDEs driven by a Brownian motion. This 
control was derived from a moment condition in the same spirit as Kolmogorov's 
criterion for Holder continuous paths. The moment condition was verified for the 
above area by the use of known stochastic integral results, though one could also 
derive it from a construction depending on the linearly interpolated Brownian 
motion. 

There are two stages to defining the integral against a geometric multiplica- 
tive functional. The first gives a functional which is almost multiplicative (see 
[|o| for definition). The second associates, uniquely, a multiplicative functional 
to the almost multiplicative functional. 

Theorem A. 2 J7^ There is a unique geometric multiplicative functional Y 
which we call the integral of the 1-form 9 against the geometric multiplicative 
functional X . We denote this by 

Yst = I 0{Xu) dX. 



Corollary A.l One has the following control on the p-variation ofY: 



e{Xu) sx 



< [Cu;{s,t)]^' /[P[z/p]\] 1 = 1,. ..,[p] (34) 



where C depends on p, ||/||Lip('y), 7, A, /3, L and [p]. 

The estimate is derived from estimating both the almost multiplicative func- 
tional and the difference of it from the integral. 

We now state two lemmas which help prove that the solutions of (|l|) are 
homcomorphic flows when the initial condition is varied. 

Lemma A.l Let X be in OG(K'^)p controlled by a regular ujq. Let f : R" 
hom(M'^,K") be a Lip(7) map for some 7 > p. Let Y}f \ i = 1,2 denote the 
element in r2G'(R")P which solves the rough integral equation 

nt 

Yst^ = / f{Y'^''>)SX 



with initial condition Yq = Oi, i = 1,2. Let Wgt be the multiplicative func- 
tional which records the difference in the multiplicative functionals Y^l^ and 
Y^^^ \ Then 

where 9 = \ai —02], co = Cujq, the constant C depends on p, ||/||Lip(7)7 /3, 7- 
The bound holds for all times s < t on the interval J = {u : w (0, u) < 1}. 
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Lemma A. 2 With the assumptions of Lemma A.l one can estimate the differ- 
ence of the increments of y}P and for any pair of times < s < t which 
satisfy lo (s, t) < 1 os follows: 



(1) 



1 



_/3(l/p)! 

In particular for any t > one has: 



(lu (0,s)+u;(0,s)(i/'')) 



LO is,t) 



(1/p) 



/3(l/p)! 



(1) 



r/'^l < \ai-a2\ C{t). 



(36) 



Now we can prove that the solutions form a flow of homeomorphisms as the 
initial condition is varied. 



Proof of Theorem 1.3 The continuity of solutions follows from Lemma A. 2 



It remains to show that the inverse map exists and is continuous. This can 
be checked by repeating all the previous arguments using the reversed path 

{Xt-s)o<s<t as the integrator. □ 

The induction part of the proof of Lemma A.l will require the following 
lemma about rescaling: 



Lemma A. 3 JT^ Let X he a multiplicative functional in r([pl)(R'*) which is of 
finite p-variation controlled by lu. Let (X, Y) be an extension of X to T'-[^1^(]R'*® 
R") of finite p-variation controlled by Klu. Then {X,(j)Y) is controlled by 



max|l,0'=''/'X ■.l<k<i<[p]juj 
where G M. In particular, if 4> < K^^^/'p < 1 then {X,(j)Y) is controlled by uj. 



Proof of Lemma 



A.l 



We set up an iteration scheme of multiplicative function- 
als which we will bound uniformly, by induction. A projection of the sequence 
proves that a Picard iteration scheme converges to the solutions of starting 
from ai and 02- Another projection shows that the difference of these solutions 
is bounded. 

Let e > and 77 > 1. Let V^l^ be the geometric multiplicative functional 
given by 

T.(l) A .7(1)(1) ^(1)(1) v(l)(0) y(2){l) ^(2)(1) ^(2)(0) 



f{ai)6X-ai,J f{ai)6X,ai,J f{a2)5X~a2, 
f{a2) SX, 02, ^ /(ai) - f{a2) 6X, e^^X.^^. 
The iteration step is a two stage process. Given we set 



Path-wise solutions and Levy processes 



26 



where is the 1-form on ((R")®'^ ® R'^) given by 
fc™(ai, ... ,08) {dAi, . . . ,dAs) 

ai g(a2,a3) dAs, dAs + j]^"'dAi, dA2, a4 g{a5,ae) dAs, 

dAe + T]~"^dA4, dA^, g{a2, 04) dAs, dAs 

g{x,y) is the 1-form appearing in Lemma 3.2] wliicli satisfies the foUowing 
relation witli respect to /; 

rix)-r{y)^Y.('^-yy g^={x,y). 



y(m+i) jg .^gjj (iegj^ed because g and fc™ are botli Lip(7) for some 7 > p — 1. 

We define l/t^+i) to be the geometric muhiphcative functional obtained 
by rescaling the first and fourth components of by erj and the seventh 

component by e. 

The uniform bound on the iterates {V^™'')m>i wiU be obtained by induction. 
X is controUed by a regular ujq so there exists a constant C such that V'^^'^ is 
controlled by w = Cluq- Suppose that V^''\k < m) are controlled by uj. From 
(Corollary A.l) there is a constant Ci such that is controlled by Ciuj. 

If we choose e > 0, > 1 such that e <Ci'^ and erj < Ci'' , then Lemma |A.3| 
implies that is controlled by uj, completing the induction step. 

The uniform control on the iterates ensures the convergence of 

|y(i)(m)}^^^^^ to the solutions of 

dr/'' = /(r/^)) d^t Y^'^=a,, z = l,2. 

Through the definition of , the sequence at the level of the paths 

will converge to the scaled difference of the two solutions 9^^{Y^^^ — y^^'). For 
s, i in J one has 

which implies that 



Proof of Lemma A. 2 We define the following set of times: 

to — and = inf{w > : w u) = 1} V j £ {1, . . . , n(s)} (37) 

where n{s) = max{j : tj < s} and tn(s)+i — s . 
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We solve the differential equation starting from s and use (^5|) to show that 

where K{s) is an upper bound on the supremum over all the possible differences 
of the paths \Ys^'' — , at time s. The bound K{s) is derived recursively by 
considering the analogous upper bound for the difference of the solutions to the 
differential equation over the time interval [ti_i,ti] given below: 

\yI'^ - I < - Yi^^ I + \Wu_, u I 
which implies that 

.;(i,_i,i,)(Vp) 



Therefore 



iliT ( 

\w%\<K{t,f n 1 

.7 = 1 



f3{l/p)\ 



j 1, . . . ,n(s) + 1. 



/3(l/p)! 



< 6''= exp 



i=i 



/3(l/p)! 



w(s,t) 



(fc/p) 



/3(fc/p)! ' 



noting that uj{tj-i, tj) = 1 and using the sub-additivity of lo we obtain 



< 6''= exp 



k 



f3{i/py. 



(w(0,s)+w(0,s)(i/P)) 



w(s,t) 



(fc/p) 



/3(fc/p)! 



By considering the above bound at the level of the paths {k ^ 1) and repeatedly 
using the triangle inequality one deduces (|36|). □ 
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